Short-range particle correlations in dilute Bose gas 
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The thermodynamics of a homogeneous dilute Bose gas with an arbitrary strong repulsion be- 
tween particles is investigated on the basis of the exact relation connecting the pair correlation 
function with the in-medium pair wave functions and occupation numbers. It is shown that the 
effective-interaction scheme which is reduced to the Bogoliubov model with the effective pairwise 
potential, is not acceptable for investigating the short-range particle correlations in a dilute strongly 
interacting Bose gas. In contrast to this scheme, our model is thermodynamically consistent and 
free of the ultraviolet divergences due to accurate treatment of the short-range boson correlations. 
The equation for the in-medium scattering amplitude is derived that makes it possible to find the 
in-medium renormalization for the pair wave functions at short boson separations. Low-density 
expansions for the main thermodynamic quantities are reinvestigated on the basis of this equation. 
Besides, the expansions are found for the interaction and kinetic energies per particle. It is demon- 
strated that for the many-boson system of the hard spheres the interaction energy is equal to zero 
for any boson density. The exact relationship between the chemical potential and in-medium pair 
wave functions is also established. 



PACS number(s): 05.30. Jp, 67.40 Db, 03.75.Fi 



I. INTRODUCTION AND BASIC EQUATIONS 

The well-known experiments with the magnetically 
trapped alkali atoms [JjJ have significantly renewed in- 
terest in the theory of the Bose-Einstein condensation 
(see, e.g., Ref. Q). In particular, it has recently been 
demonstrated by the present authors Q that the cus- 
tomary way of investigating a dilute Bose gas with a pure 
repulsive and arbitrary strong interaction Q is thermo- 
dynamically inconsistent. At n = N/V — > this way is 
known || to be reduced to the Bogoliubov model Q with 
the "bare" pairwise potential <I>(r) replaced by an effec- 
tive, "dressed" one. This is why below the approach of 
Ref. Q is called the "effective-interaction method" . The 
"dressed" pairwise potential is usually derived by sum- 
ming the ladder diagrams and involves, as is assumed, 
all the necessary information on the short-range spatial 
correlations of bosons Q. In the final expressions use 
of the effective interaction results in substituting the ex- 
act scattering length a for its Born approximation ao [pi . 
This allows for operating with strongly singular poten- 
tials, but at the price of loss of the thermodynamic con- 
sistency. On the contrary, the strong-coupling general- 
ization of the Bogoliubov model proposed by the present 
authors in the paper || , is based on the variational pro- 
cedure and does not invoke any mean- field arguments. 
Owing to this structure of the generalization we do not 
need to worry about the thermodynamic consistency. 

The trouble mentioned above gives rise to various 
misrepresentations of the effective-interaction approach. 
For example, the condition of the self-consistency leads 
to zero condensate depletion within the pseudopotential 
model j?j . Another its manifestation is an irrelevant pic- 
ture of the pair boson correlation at short particle sepa- 
rations. This important point calls for a comprehensive 



analysis which has not been fulfilled in Ref. || for reason 
of space. Thus, in the present paper we continue rein- 
vestigation of a dilute Bose gas with an arbitrary strong 
repulsion between particles within the model proposed 
in the paper ]3| , the short-range boson correlations being 
of special interest now. The zero temperature is under 
consideration below. 

The formalism of the present paper is concerned with 
the reduced density matrix of the second order (the 2- 
matrix) and its eigenfunctions which we call, following 
Bogoliubov [gj, the in- medium pair wave functions. As 
the 2-matrix with its eigenfunctions are not often dis- 
cussed in the modern scientific literature on the Bose- 
Einstein condensation, it is worth noting some basic no- 
tations and formulas. The 2-matrix for the many-body 
system of spinless bosons can be represented as || : 



p 2 (ri,r^;ri,r 2 ) 



F 2 (ri,r 2 ;ri,r 2 ) 



N(N-1) 

where the pair correlation function is given by 
^ 2 (n,r 2 ;ri,r 2 ) = (^t( ri )^t( r2 )^( r ^( r ^). 



(1) 



Here ip(r) and ip^r) denote the boson field opera- 
tors. Use of the pair correlation function that dif- 
fers from the 2-matrix only by the normalization fac- 
tor, is more convenient in the thermodynamic limit (n — 
N/V — const, V — ► oo) when F 2 (r l5 r 2 ; r[, r' 2 ) oc I 
while p 2 (r^ , r,; ri, r 2 ) oc 1/V 2 . Recently it has been 
found |3 10 Tlj] that for the uniform Bose gas, the sys- 
tem with a small depletion of the zero-momentum state, 
the correlation function ([!]) can be written in the ther- 
modynamic limit as follows: 

F 2 (r 1 ,r 2 ;ri,r / a )=^*(rMr') 



1 



/d 

x exp[iq- (R'-R)], (2) 

where r = ri — R = (ri +T2)/2 and similar rela- 
tions take place for r' and R', respectively. In Eq. (g) 
no = Nq/V is the density of the particles in the zero- 
momentum one-boson state, n q — (a q a q ) stands for the 
distribution of the uncondensed bosons over momenta, 
<^(r) and <£> q /2(r) are pair wave functions in medium in 
the centre-of-mass system. Namely, <p(r) is the wave 
function of a pair of particles being both condensed. In 
turn, </? q / 2 (r) denotes the wave function of the relative 
motion of a pair of bosons with the total momentum ?iq, 
this pair including one condensed and one uncondensed 
particles. So, Eq. (0) takes into account the condensate- 
condensate and supracondensate-condensate pair states 
and is related to the situation of a small depletion of 
the zero-momentum one-boson state. For the pair wave 
functions we have 

p(r) = l+V(r), Mr) = \/2cos(p-r) + ^ p (r)(^0), ( 3 ) 

where the scattering waves ip(r) and ipp(r) obey the 
boundary conditions for r — > oo: 



-► 0, %{r) -» 



(4) 



which follow from the Bogoliubov principle of correlation 
weakening The Fourier transforms of the functions 
tp(r) and ipp( r ) can explicitly be expressed in terms of 
the Bose operators at and a p pTt : 



ip(k) = (akS-k)/no, ^p( k ) 



V (a 2p a p+k a p _ k ) 



2n 



n 2p 



(5) 



In the representation (^|) the terms corresponding to 
the supracondensate-supracondensate "channel" , are ne- 
glected, i.e. we omit the contribution of pairs of the 
particles which are both uncondensed. Besides, it is as- 
sumed that there are no bound states of pairs of bosons, 
which is obviously realized for a purely repulsive interac- 
tion between bosons. 

The diagonal matrix element of the pair correlation 
function (0) is proportional to the pair distribution func- 
tion 



g(r) = F 2 (ri,r2;ri,r 2 )/n 2 



(6) 



that can directly be observed in the scattering experi- 
ments. Derivation of Eqs. (||)-(||) and detailed discussions 
can be found in Ref. jl(J . 

The two limiting cases n — > and r — > correspond 
to the situation when the behaviour of two particles in 
medium is determined by the ordinary two-body prob- 
lem provided the pairwise interaction <&(r) is repulsive 
and goes to infinity at short boson separations. In par- 
ticular, when n — > we have (n — no)/n — > and, as it 



is known since the Bogoliubov original paper and follows 
also from Eqs. (|^) and (^|), 



ff(r)-b<°>(r)] 



(7) 



Here ip(°\r) is defined by ip^\r) = limn—,0 and 
obeys the ordinary two-body Schrodinger equation in 
the centre-of-mass system (g) (see Ref. ||). Similarly, 
(fip^ (r) = lim n ^o ( r ) obeys the Schrodinger equation 
related to the eigenvalue h 2 p 2 /m that corresponds to the 
relative motion of two particles with the momentum hp. 
These conditions should be satisfied in any theory that 
appropriately takes into account the short-range corre- 
lations of particles. Below we show that our model || 
leads to the correct picture of the spatial correlations in 
contrast to the effective-interaction approach that leads 
to negative values of the pair distribution function (see 
Sec.0). 

Having at our disposal the distribution function n^ and 
the set of the pair wave functions ip (r) and </2 p (r), we are 
able to calculate the main thermodynamic quantities of 
the system of interest. In particular, the mean energy 
per particle is expressed in terms of n^ and g(r) via the 
well-known formula (see, e.g. Ref. M) 



(8) 



where T fc = h 2 k 2 /2 m is the one-particle kinetic energy, 
m is the "bare" mass of particles, and n = N/V stands 
for the boson density. 

The organization of this paper is as follows. In Sec. |j| 
we give, for convenience, helpful information concern- 
ing the classification of the pairwise potentials used in 
the ordinary two-body problem. In Sec. Ill the Bogoli- 
ubov model of a weakly interacting Bose gas is consid- 
ered within a variational scheme. This scheme yields the 
system of two equations connecting n^ with <p(r). As 
to the supracondensate-condensate pair wave functions, 
they are the symmetrized plane waves in the Bogoliubov 
model: ip p (r) — 0. In the next section the low-density ex- 
pansions for the condensate depletion and mean energy 
per particle of a weak-coupling Bose gas are calculated 
within the Bogoliubov model. The effective-interaction 
approach of Ref. j| is analyzed in section [v|. Using the 
results of the previous sections [H and |iv| , we show that 
the effective-interaction approach is thermodynamically 
inconsistent. This inconsistency turned out to be di- 
rectly related to an irrelevant picture of the short-range 
spatial boson correlations. In particular, for a strongly 
singular potential the effective-interaction scheme yields 
for the pair distribution function the nonphysical result 
g(r = 0) = — 1 in the limit n — > 0. It is also demonstrated 
that the well-known ultraviolet divergence appearing in 
the effective-interaction approach as well as the ther- 
modynamic inconsistency occur because the Bogoliubov 
framework is used beyond the range of its validity. The 
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regularizing procedure which consists in omitting the 
divergent integral J d 3 fc/fc 2 , can be justified provided 
the quantities of interest depend on the pairwise poten- 
tial through the mediation of the scattering length (pj]). 
Sec. VI regards a correct strong-coupling generalization 
of the Bogoliubov model. This generalization is based on 
Eq. (H) taken together with Eq. (g) rather than on its lin- 
earized variant Eq. (|25| ) used in the effective-interaction 
approach and being the weak-coupling approximation for 
g(r). A variational procedure similar to that of Sec. [II 
is formulated. It provides the system of equations which 
should be solved to find the pair wave functions in con- 
junction with the momentum distribution. For a dilute 
Bose gas this system is reduced to the set of two equa- 
tions connecting rik and ip(r). There is essential differ- 
ence between these equations and those of Sec. III. Now 



the pairwise potential $(r) appears only in the combi- 
nation tp(r)$(r) which allows for using the strongly sin- 
gular potentials beyond the effective-interaction scheme. 
In what concerns the supracondensate-condensate contri- 
bution to the thermodynamic quantities, it can be calcu- 
lated with the relation lim p ^o¥>p( r ) — V% ( p(j') resulting 
from Eqs. (|j) and In the next section we investigate 
the short-range renormalization for tp(r) conditioned by 
presence of surrounding bosons. Its long-range behavior 
is also discussed. All the investigation of this section is 
based on the in-medium Lippmann-Schwinger equation 
coming from the equa tion for ip(r) found in the previous 
section. In Sec. VIII the low-density expansions for the 
Bose condensate depletion, the energy per particle and 
the chemical potential are found within the model pre- 
sented in Sec. various calculation ways being used. 
For this purpose we establish the exact relationship be- 
tween the chemical potential and pair wave functions in 
a condensed many-boson system. Here we also evaluate 
the kinetic and interaction energies per particle which, 
to our knowledge, have never been calculated. It should 
be stressed that they explicitly depend on a shape of the 
pairwise potential even in the leading order of the low- 
density expansion. In the framework of our approach we 
are able to perform all the calculations concerning the 
kinetic and potential energies both directly and with the 
Hellmann-Feynman theorem, in contrast to the effective- 
interaction method. Main results and prospects are dis- 
cussed in the last section. 



II. CLASSIFICATION OF INTERACTION 
POTENTIALS 



Before further consideration we recall the classification 
of the pairwise interactions $(r) which is used in the ordi- 
nary two-body problem. In this paper we only deal with 
the short-range potentials which go to zero for r — > oo as 
$(r) — > l/r m (m > 3), or even faster. Let us consider 
the solution of the two-body Schrodinger equation in the 
ccntre-of-mass system 



-VV°>(r) + *(r)¥>W(r) = 



(9) 



which corresponds to the scattering state with the mo- 
mentum p — 0: ip(°'(r) = 1 + ip(°\r), where the scatter- 
ing part behaves as 



V/ 0) (r) _ a / r 



(10) 



when r — > oo. Owing to this boundary condition with 
the real quantity a, the solution y>'°)(r) is chosen to be 
real also. The scattering length a is defined by means of 
the scattering amplitude EA°)(0): 



rf/ (0) (0), 



(11) 
(12) 



As applied to Eq. (|9|), the perturbation technique gives 
the expansion for its solution 



Airh 2 

[/(°)(0) = / d 3 rV )(r)$(r). 



0(°>(fc) = ^ 0) (fc)+7/4 0) (fc) + 



(13) 
(14) 



which leads to the following expansion for the scattering 
length (O): 



a = a a + ai + a,2 + 
rn 



a Q = 



:<f>(k = 0), oi 



d 3 fc $ 2 (fc) 
Airh 2 J (2tt)3 2T k ' 



(15) 
(16) 



Here ijj^°\k) and $(fc) stand for the Fourier transforms 
of ip(°\r) and ^(r), respectively. If we restrict ourselves 
to the first terms in Eqs. © and @ h/> (0) (fc) ~ ^ { °\k) 
and a ~ ao] we arrive at the Born approximation for the 
wave function and the scattering length, respectively. 

The interaction is called the weak-coupling one pro- 
vided the Born approximation works well, in particular, 



d 3 fc $ 2 (fc) 
(2tt)3 2T fc 



<$(fc = 0). 



(17) 



This is valid if, first, the potential is integrable, and, 
second, it is proportional to a small parameter, the cou- 
pling constant. The latter implies that |V^°H r )l "C 1, 
and, so, the Born approximation (|l4|) is nothing else but 
linearization of Eq. (||) with respect to ij}^(r) propor- 
tional to the coupling constant: 



V 2 (l+^°)(r)) 
(l + #°)(r)) 



(0) 



(r). 



The potential is called singular if it is integrable but the 
Born approximation does not work well. At last, the 
potential is of strongly singular, or hard-core, type if it 
is not integrable [$(r) — > l/r m (m > 3) for r — ► 0], 
and the terms ([lj) and © ) cannot thus exist. In the 
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present paper the pairwise interaction of this type is ex- 
actly implied when we speak about the strong-coupling 
regime. For example, the well-known Lennard-Jones po- 
tential corresponds to this case together with the hard- 
sphere interaction 



$(r) 



0, 



r < a, 
r > a. 



(18) 



In the strong-coupling regime the solution of Eq. (0) 
obeys the boundary condition (p^(r = 0) = 0, otherwise 
the interaction energy £ m( = j d 3 r {(p^ (r)] 2 $(r) and 
the scattering length (|ll]) would be infinitely large. 

In further consideration we make use of the variational 
theorem for the scattering amplitude fll^): 



«/ (0) (0) 



d 3 r 



^(0) (r)(5 (_!Lv 2 W )(r) 

V TO / 



(19) 



In order to prove this relation, we represent Eq. ( |12| ) in 
the form 

r vfi 2 i 
C/ (0) (0) = / d 3 r — |V^ (0) W| 2 + [v {0 Hr)] 2 $(r) , (20) 
J Lto " J 

which can be found using integration by parts and tak- 
ing into account the Schrodinger equation (Q) and the 
boundary condition (|l0). Further, varying Eq. ( pp[ ) and 
keeping in mind Eqs. (^) and ([l0|), we arrive at Eq. (19|). 
The relation (19) bears analogy to the variational theo- 
rem for the energy: 



SE n 



d 3 r^°)(r)j(--V 2 + $(r))^°)(r), 

V TO / 



here the real wave function ipn\v) obeys the Schrodin- 
ger equation for a bound state with the energy E n < 0. 
Equation (|2^) can be represented in the more convenient 
form 



d 3 r [v3 (0) (r)] 2 $(r) = Airh 2 {a - b)/m. 



(21) 



Here one more characteristic length b (in addition to a) 
has been introduced: 



1 

4-7T 



d'V |VV> (0) (r) 



(22) 



It follows from Eq. (£^) that b is a positive quantity. We 
stress that b is not expressed in terms of a and depends on 
a particular shape of the interaction potential <&(r). For 
example, when $(r) is the hard-sphere potential (Jig), we 
have b — a. While for <I>(r) close to zero, in the weak- 
coupling regime, we have b ~ — ffli, a — ao j and, hence, 
b < a. 

Lastly, from the definitions ([ll]) and ( p2| ) and the vari- 
ational theorem ( [l9|) it follows that 



7 



da 



da 
dm 



(23) 



where we introduce the auxiliary parameter 7 called cou- 
pling constant [i.e. <I>(r) — > 7$(r)]. The first equality in 
Eq. (^3| ) implies that a = a('ym), which is an obvious 
consequence of the definition (|il| ) and the Schrodinger 
equation (5h. The relations fl23|) demonstrate that the 
quantity b is expressed in terms of a and its derivative 
with respect to 7 (or m) rather than in terms of a. 



III. VARIATIONAL TREATMENT OF THE 
BOGOLIUBOV MODEL 

Although the aim of this paper is to investigate a di- 
lute Bose gas with the strong-coupling interaction, it is 
deductive to start with the Bogoliubov model related to 
the weak-coupling regime. This regime implies a minor 
role of particle scattering, both in medium and out of it, 
and, thus, is characterized by the following inequalities 
for the scattering waves (J3J) : 



|V(r)| « 1, |^p(r)| « 1. 



(24) 



In particular, the Bogoliubov model operates with the 
choice 

|V(r)|«l, V P (r)=0. 

As the depletion of the Bose condensate (n — no)/n 
is small in a weakly interacting many-boson system, 
we have for the one-particle density matrix F\ (r) = 
(^(r^Cra)): 



F 1 (r) 




J 


n 







d k nk 
(2?r) 3 ~n 



exp(ik • 



< 



n 



< 1. 



So, the Bogoliubov scheme of treating a Bose gas involves 
two small quantities ip(r) and Fi(r)/n, and completely 
neglects scattering in the supracondensate-condensate 
sector of g(r): ip p {r) — 0. This along with Eq. (||) al- 
lows for rewriting Eq. (||) in the following form: 



g(r) = 1 + 2ip(r) + - 
n 



d 3 fc 
(27p 



nfcexp(ik-r). (25) 



Here we have restricted ourselves to the terms linear in 
V>(0 and Fi(r)/n. Besides, it is implied that ip*(r) — 
ipir), for the pair wave functions can be chosen as real 
quantities. Inserting Eq. (|25|) into Eq. (||), we are able 
to employ a variational procedure to derive the unknown 
quantities ip{k) and n^. In so doing, we should realize 
that nk and tp(k) are not independent variables. In fact, 
there are no spatial boson correlations in absence of inter- 
action Hence, in this case ^>{k) = 0, and, as we in- 
vestigate the ground state, all the bosons are condensed, 
iik = 0. While in presence of interaction ip(k) ^ 0, which 
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leads to a nonzero depletion and, so, n k ^ 0. Within the 
Bogoliubov model ip(k) is related to n k by 



n k (n k + 1) = nlip 2 (k). 



(26) 



Indeed, according to the canonical Bogoliubov transfor- 
mation, quasiparticle (bogolon) operators d^ and dk are 
connected with the operators of the primordial bosons by 
the expression 



a k = tifed k + Vfcdt k , 
where 

ui - vi 



T-i = Uk&' k + v fe d_ k , (27) 



1. (28) 



Within the Bogoliubov model the ground state of the 
system of interest is the bogolon vacuum, and, so, at the 
zero temperature we have 

(4<*k) = °- ( 29 ) 
Then, using (g), @ and (f|), one can find 

n k = vl, tp(k) = u k v k /n , 

which in conjunction with ( psj ) leads to Eq. (p6|) . We re- 
mark that beyond the Bogoliubov scheme Eq. (|26|) is not 
valid and should be corrected (see Eq. ( p5l| ) in the present 
paper and discussion on this question in Ref. ||). 

Now, inserting Eq. (|25|) into Eq. (@) and varying the 
obtained expression with respect to ip{k) and n k , we de- 
rive 



- 2T k i\)(k) = $(fc)(l + 2n k ) + 2r#(fc)*(*). (33) 
which, with the help of Eq. (p5|), can be represented as 

-VV(r)=*(r)+n / d 3 y *(»)(«/(|r - y|) - l). (34) 

TO J 

Equation (|34|) is very similar to the Bethe-Goldstone one. 
Necessary details concerning Eq. ([54]) can be found in 
the papers PLR- We remark that the right-hand side 
(r.h.s.) of Eq. fl34|) can be thought of as the in- medium 
potential of the boson-boson interaction in the Born ap- 
proximation. Indeed, Eq. ( [34]) is derived from the more 
general equation given by Eq. ( |56| ) by means of lineariza- 
tion in ip(r) and F\(r)/n. So, the Bogoliubov model can 
be treated as the in-medium Born approximation for its 
generalization developed on the basis of Eqs. @\ and (||) 
beyond the weak-coupling regime (see section fVU) . In ac- 
cordance with this treatment, Eq. (B4) at n = is nothing 
but the Fourier transform of Eq. (|14|) while Eq. (|5^) is 
reduced to the exact Schrodinger equation (g) at n = 0. 
We will return to this important point in Sec. [v| 

The system of Eqs. ( p6| ) (here we should set n = no) 
and ( p3h can easily be solved, which leads to the familiar 
results ifll : 



1 



n k 



ip(k) 



T k + n<S>(k) 



2 \ y/T* + 2nT k $(k) 

1 $(fc) 

2 y/T£ + 2nT k ${k) ' 



- 1 



(35) 
(36) 



Se 



d J fc r 



(2tt) 3 



(r k + n*(Jfe)) — + n$(k)6i/>(k) . (30) 



According to Eq. ( |26| ) infinitesimal changes Sip(k) and 
5n k are connected by 



(2n k + l)Sn k ip(k) f d 3 q 



2nli>{k) 
where the equality 



no 



no 



W) 



(27T) 



3 5n q , 



(31) 



?n k 



is implied. Taking 5e — and using Eqs. (30) and (|3 
we find the following equation: 



2T k i/j(k) = -5$(fc)(l + 2n k ) + 2m/j(k) 



x $(Jfe) 



n 

n 



d J q 



*(#(?) • (32) 



Note that Eq. ( p2[ ) is able to yield results being accurate 
only to the leading order in (n — no)/n because Eq. (25) 
is valid to the next-to- leading order So, Eq. (32) 

should be written in the form 



IV. DENSITY EXPANSIONS IN THE 
BOGOLIUBOV MODEL 

As it was mentioned, in this paper we investigate the 
strong-coupling regime for a dilute Bose gas. So, con- 
sidering a dilute Bose gas in the weak-coupling approx- 
imation can be a good exercise providing us with useful 
information. Let us investigate the thermodynamics of a 
dilute many-boson system within the Bogoliubov model. 
With Eqs. (§), @, @ and @ we derive 



'Mo) 



i 

2n 



d*q 

WT 3 



x (y/l* + 2nT q <P(q) —T q — n$(q)) . (37) 



Here we describe in details the method for obtaining the 
low-density expansions for expressions like Eq. (|37| ) . This 
equation can be represented in the following form: 



e = - ( $(0) - 



d 3 g $ 2 (g) 
(2^) 3 2T q 



(38) 



where 



5 



I = 



A 6 q 
(2+j 3 



T' 2 



T T 
n n 



<& 2 (g) 

2(T,/n) 



Now, with the "scaling" substitution 
q = q'v2mn/?i 



we derive 



,3/2 



2m 



3/2 



3 „' 



:f(q\n) 



(39) 



(40) 



where 



/(<?', ™) = yV) 4 + 2(<7') WV2W?i) - (g') 2 



-$(q'V2mn/fr) + 



2(9') 



A2 



The advantage of the representation ( |38| ) is that the 
resulted integral in Eq. ( p0| ) uniformly converges for 
q' — > oo with respect to n for n — + 0, and, thus, we 
obtain 



lira 



o n 3 / 2 



2m 



3/2 



dV 

(2+1 3 



/(g> = 0). 



Here the integral is readily calculated, and the main 
asymptotics for I(n) is given by 
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Further, with the help of Eqs. @, @ and @, the 
expression (|37|) is rewritten as 



2irh 2 n(ao + a±) 2i:h 2 nao 128 



15^ 



(42) 



Thus, we obtain the first two terms in the density expan- 
sion for the energy per particle within the Bogoliubov 
model. 

The density expansion for the condensate depletion is 
inferred from Eq. (35) in the same manner employing the 
substitution 



c 



n — hq 



n 



J 3 

d q n q 
(2tt) 3 n 




(43) 



Now let us discuss the range of validity of the Bogoli- 
ubov model. First, the condensate depletion ( |43"| ) should 
be small as the representation for the pair correlation 
function which we start from is valid only in this 
case. Note that if we expand the depletion with respect 
to the coupling constant 7 (we assume that $(r) cx 7), we 
also arrive at Eq. ( |43| ) since the occupation number ( |35| ) 
depends only on the production n~f. Thus, the condition 



na,Q <C 1 should be fulfilled. Second, we exploit the weak- 
coupling character of the pairwise interaction $(r), which 
implies that the condition (17) should also be satisfied. 
Note that Bogoliubov himself realized this necessary con- 
dition since he treated the term n$(0)/2 = 2nh 2 nao/m 
involved in the mean energy ( |37j ) as the major one |J. 
Beyond the inequality ( |l7| ) the model may be thermody- 
namically unstable. In particular, the opposite case 



$(0) < 



d 3 fc $ 2 (fc) 
(2tt) 3 2T k 



(44) 



leads to the negative scattering length in the next-to- 
Born approximation a = ao + &i < which at sufficiently 
low densities results in the incorrect sign for the com- 
pressibility —d 2 E/dV 2 = dP/dV > as it is seen from 
Eq. (^||). Note that this important point is not always 
stressed in literature. Moreover, Bru and Zagrebnov |15[| 
proposed the model reduced to the Bogoliubov approach 
but in conjunction with the inequality (Q). We have to 
conclude that this model hardly has a physical sense. 



V. SHORT-RANGE CORRELATIONS AND 
ULTRAVIOLET DIVERGENCE WITHIN THE 
EFFECTIVE-INTERACTION APPROACH 

After the detailed consideration fulfilled in the previous 
sections, we can argue that the Bogoliubov model with 
the "bare" potential <&(r) replaced by an effective one 
Q, is thermodynamically inconsistent. Indeed, the basic 
relations of the Bogoliubov model (25) and (|2g) depend 
on interaction implicitly. Hence, the pairwise potential 
appearing in Eqs. (|3^)-(|3^) comes from Eq. (|) that is 
the general relation valid in both the weak- and strong- 
coupling regimes. Thus, a calculating procedure based on 
@ and @ has to eventuate in Eqs. Other- 
wise, like in the case of using the Bogoliubov model with 
an effective interaction, this procedure does not yield the 
result minimizing the mean energy. Note that we do not 
mean, of course, that the t-matrix approach or the pseu- 
dopotential method can not be applied in the quantum 
scattering problem. It is only stated that the usual way 
of combining the ladder diagrams with the RPA ones 
( "bubbles" ) leads to the thermodynamic inconsistency. 

To clarify the reason for this inconsistency, let us take 
a look at the picture of the spatial boson correlations 
derived in the framework of the effective-interaction ap- 
proach. According to the paper by Hugenholtz and 
Pines [Q (see Eq. (5.10a) therein), the structural fac- 
tor 



S(k) =l + n / d 3 r [g(r) - 1] exp(-ik • r) 



(45) 



of a strong-coupling Bose gas can be written at n — > as 
follows: 



S(k) = 1 + 2^(4a k ) + ^«a k a- k ) + (a! k a k )) . (46) 



no 
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Using Eq. the equality ip(k) = tp*(k) Jl6| and defini- 
tion (^5|) of the structural factor, one can readily verify 
that Eq. @ is reduced to Eq. @. This Bogoliubov 
relation does not depend on the interaction potential ex- 
plicitly. So, use of the "dressed" interaction can in no 
way disturb the form of Eq. (|25|), and, therefore, the 
effective-interaction approach deals with the pair distri- 
bution function whose structure has the obvious weak- 
coupling character. In particular, from Eq. ( |46| ) it can 
be found that g(r) — > 1 + 2i)j^\r) for n — > as opposed 
to the correct strong-coupling result given by Eq. (|7]). 
However, the wave function (r) obeys Eq. @ in the 
effective- interaction approach JL7J . While within the Bo- 
goliubov model (p(°)(r) is the solution of Eq. (Q) at 
n = 0. This equation (Eq. (Q) at n = 0) comes from the 
Schrodinger equation (|^) in the Born approximation (see 
discussion in Sec. lll|). Thus, the effective- interaction ap- 
proach is not totally reduced to the weak-coupling frame- 
work due to its features of the strong-coupling character. 
Exactly the combination of the peculiarities of both the 
strong- and weak- coupling regimes is the reason for the 
thermodynamic inconsistency mentioned above. 

It is also worth noting that this combination of the 
features of weakly and strongly interacting Bose gases 
leads not only to the thermodynamic inconsistency. It 
results also in one more misrepresentation of the effective- 
interaction approach. We mean an irrelevant picture of 
the short-range boson correlations. Indeed, in the case 
of a strongly singular pair interaction for the solution of 
Eq. (Q) we have <p(°J(r = 0) = [see Sec. [ll], which pro- 
vides -0^°' (r = 0) = — 1. Within the effective-interaction 
scheme g(r) obeys Eq. (p5[) while t^ ' (r) satisfies Eq. (^). 
This implies that g(r = 0) -> 1 + 2ip(°>(r = 0) = -1 in 
the zero density limit when (n—no)/n — > 0. The obtained 
result does not agree with the physical sense of g(r) (the 
conditional probability) and has nothing to do with the 
strong-coupling regime when the relation g{r = 0) = 
has to be satisfied. The situation even aggravates if we 
recall that the scattering parts of the supracondensate- 
condensate pair wave functions tpp( r ) are equal to zero 
in the Bogoliubov model. So, in what concerns the pair 
distribution function, the "triple" correlations involved 
in Eq. (|^) for ipp(k) are completely ignored within the 
effective-interaction scheme. However, when deriving an 
equation for the effective potential, these correlations are 
taken into consideration, for example, within the Beliaev 
approach (see discussion in the review jl8|). So, we face 
one more combination of the weak- and strong-coupling 
features being characteristic of the approach of Ref . jl| . 

Since the short-range behaviour of the pair distribution 
function is not correct within the effective- interaction ap- 
proach, one can expect some problems related to evalua- 
tion of the mean energy (||). Let us consider the effective- 
interaction method in its simplest variant, the so-called 
pseudopotential model (see paper by Lee, Huang and 
Yang of Ref. JiJ^. This variant implies the replacement 
$(r) — > <5(r)47r?i a/m, and, hence, for the Fourier trans- 
form we have 



$(fc) — > 4irh 2 a/m = const, 



(47) 



where a is the scattering length ( |ll| ) obtained from the 
Schrodinger equation (^). So, the pseudopotential model 
is reduced to the Bogoliubov one with the effective pair- 
wise interaction given by (|47|). In the well-known text- 
book g one can find two ways of calculating the leading 
and next-to-leading terms of the low-density expansion of 
the energy of a dilute Bose gas within the pseudopoten- 
tial model. One of them, see the pages 314-319, consists 
in dealing directly with the Hamiltonian of the system 
and faces the divergent integral J d 3 k/k 2 (the ultravio- 
let divergence). The latter, given on the pages 218-223, 
allows for calculating the difference e — fi/2 and does not 
lead to any divergence. 

In the previous section we have derived the low-density 
expansion ( p2[ ) corresponding to the Bogoliubov model. 
This expansion can help us to understand reasons for the 
ambiguous result of the pseudopotential model. Use of 
the pseudopotential ( |47| ) leads to the substitution a — > ao 
in Eq. (|42j). Besides, a\ — > —00 as it becomes propor- 
tional toj d 3 k/k 2 [see Eqs. (|l6|)]. This agrees with the 
result of the first way of calculating e in the textbook @] . 
The divergent integral is usually removed because it is 
assumed that "this divergence is not very basic" ||. So, 
we arrive at the correct expression ( |70| ) which is found 
in our model beyond any divergences. The reason for 
the singularity is obvious because the necessary condi- 
tion of the validity of the Bogoliubov model is not 
satisfied. However, the question remains why the pseu- 
dopotential approach results nevertheless in the correct 
final expression ([70|)? The point is that the effective- 
interaction scheme actually involves an additional as- 
sumption, namely, the Landau postulate (see the foot- 
note in Ref. || and discussion in Ref. j|]). This postu- 
late asserts that the properties of dilute quantum gases 
are ruled by the scattering length a |19|. Let us consider 
how the additional assumption is used when deriving the 
low-density expansion for the mean energy. According 
to the Landau argument this expansion should be of the 
form 



d(a)n + c 2 (a)n 3/2 + 



(48) 



where the factors Ci can depend on various quantities but 
one of them, the "bare" potential $>(r), is involved only 
through the mediation of the scattering length (|l|) . Sub- 
stituting the Born series (|l5|) in the expressions for c, (a) 
in the weak- coupling regime, we obtain 



Ci(a) = c,;(a + a\ H ) 

dci(a) 



Ci(a ) + 



do 



a 1 



(49) 



As the functional dependencies Ci(a) are of the same form 
in both the weak- and strong-coupling regimes, one is 
able to restore them by keeping the Born terms Cj(ao) 
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in the expansion (Q) and omitting others (dependent 
on Oi, 02, •••). It can readily be verified in this way 
that Eq. @ leads to Eq. ©. Thus, the pseudopo- 
tential approach provides (after the regularization) the 
correct result given by Eq. ([f0|) because it is equivalent, 
in the first two orders of the low-density expansion, to 
the calculating scheme using the Bogoliubov model to- 
gether with Eqs. (^8|) and ( f49|) based on the Landau 
postulate Note that this simple scheme looks even 
more accurate and justified than the pseudopotential ap- 
proach. At least, it allows for investigating a strongly 
interacting Bose gas beyond any ultraviolet divergence 
which appears as a result of violating subtle balance of 
the correlation terms coming from the boson-boson scat- 
tering. However, both the pseudopotential approach and 
the Bogoliubov model used together with Eqs. (48) and 
( [49| ) cannot yield adequate microscopical results concern- 
ing the strong-coupling regime. 

The second way of calculating e(n) within the pseu- 
dopotential model allows one to find the low-density ex- 
pansion (Ptq) starting from the difference 



= e — 



1 d(en) 

2 dn 



2-kTi nao 32 



(50) 



where Eq. and the well-known thermodynamic rela- 
tion a = d(ne(n)) /dn are of use. There is no divergent 
integral here due to the specific property of the expan- 
sion (f42f): ax is only involved in the leading-order term 
being exactly cancelled in Eq. ( |50| ) p5| . The solution of 
the differential equation (pO) (after replacing ciq by a) 



is of the form e 




To spec- 



ify ci being a constant of integration, one again needs 
to involve information additional to Eq. (|50|). Following 
Landau and postulating that ci depends on the pairwise 
potential only through the mediation of the scattering 
length a, one arrives at c\ — 2nh 2 a/m, which eventuates 
in Eq. @. 

Thus, we remark one more that the effective- 
interaction approach taken in conjunction with the Lan- 
dau postulate yields the correct expansion (f7(i[). Even the 
Wu's term |2l]] in the low-density expansion of the en- 
ergy of a strong-coupling Bose gas is likely to be correct 
because it is present in the weak-coupling calculations 
beyond the Bogoliubov model |22| . However, the mi- 
croscopical results found within the effective-interaction 
approach should be reexamined. So, a correct strong- 
coupling generalization of the Bogoliubov model should 
be constructed. It is also of importance that the density 
expansions for the quantities depending on the form of 
$(r), for example, the interaction ( f76[) and kinetic ( |77| ) 
energies, cannot directly be derived within the ef fective- 
interaction scheme. We discuss this point in Sec. VIII E . 



VI. STRONG-COUPLING GENERALIZATION OF 
THE BOGOLIUBOV MODEL 

To avoid the serious problems mentioned in the previ- 
ous section, we should abondon the effective-interaction 
method and use a way based on Eq. (|2|). Equations (|^), 
@ and (||) make it possible to express e in terms of the 
pair wave functions and momentum distribution. So, the 
variational procedure similar to that of Sec. Ill can be 
employed to determine these basic quantities. In so doing 
we should again keep in mind that the momentum dis- 
tribution depends on the scattering waves (see Sec. |lj). 
However, now we are not able to use the Bogoliubov re- 
lation (|2^) which does not take into account scattering in 
the supracondensate-condensate sector. In the paper || 
the following extension of Eq. (Eq) has been proposed: 



n k (n k + 1) = nlijj 2 (k) + 2n a 



(2tt) 



f n,V£/ 3 (k). (51) 



This expression has been derived with the help of the rea- 
sonable expectation that the equation for V'p(k) should 
be reduced to the equation for ip{k) in the limit p — > 0. 
It is interesting to note the obvious structural parallels 
between Eqs. (||) and (|5"l|). Now, inserting Eqs. (||) and 
(^|) into Eq. (|J) and, then, perturbing ip(k) and n k under 
the condition (pl|), we find: 

- 2T k ^{k) = U(k)(l + 2n k ) + 2ntp(k)U'(k). (52) 

Here U(k) and U'(k) are defined by 



U(k) = J d 3 r tp(r)$(r) exp(-ik • r) 



(53) 



U'(k)= / d\(vl /2 {v)-^{r))<S>{r) 



d 3 g U(q) 

(2tt) 3 ^(g) 



(^(q)-^))- (54) 



An equation for ippi^) can be derived in the same man- 
ner. Note that we have the following limiting relation: 



lim </? p (r) = V2tp(r), 



(55) 



where the factor \[2 comes from the second expression 
in Eq. (||). Using the equation for ^ p (k), one can be 
convinced that V'o(fc) = lim p _o V>p(k) = const • ip(k). We 
should put this constant equal to \/2 in order to obtain 
Eq. @ || . So, we have tp(r) - V2(p P (r) oc p 2 (see Ref. 
10 in the paper Q), which provides U'(k) — U(k) oc k A 
for k — > 0. Besides, it is easy to verify that [U(k) — 
U'(k)]/T k — * when k — * oo. Therefore at small densities 
n[U'\k) - U{k)] < T k for all momenta. This is why the 
difference between U'(k) and U(k) does not play a role 
when calculating the first two terms of the low-density ex- 
pansions for the basic thermodynamic quantities. Thus, 
at sufficiently small densities Eq. (|52|) can be rewritten 
in the following form: 
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m 

+n J d 3 y V (y)<t>(y)(g tr (\r - y|) - l), (56) 

where gt r { r ) stands for the truncated pair distribution 
function which is equal to the right-hand side of Eq. ( |25| ) 
even beyond the weak-coupling regime. 

Equations (||) [and, hence, Eqs. @, (||)] and ( |5l| ) are 
written with the condensate-condensate and supracon- 
densate-condensate "channels" taken into account. As 
these relations are accurate to the next-to-leading order 
in (n — no)/n, then, Eq. ( |52| ) can be accurate only to the 
leading order in {n — no)/n. So, it would be wrong to 
solve Eq. (jo"^) together with Eq. (|5l|). One should inves- 
tigate Eq. (p2\ ) in conjunction with the shortened version 
of Eq. ( |5l| ) given by Eq. ( |26| ) where the equality n = rio 
is implied. The system of Eqs. (^6|) and (^2|) has the 
following solution: 




m = - T 



T k + nU(k) 
f I + 2nf k U(k) 

U(k) 
f I + 2nf k U{k) 



(57) 
(58) 



interacting Bose gas of the quasiparticles with the renor- 
malized kinetic energy T k and the effective interaction 
U(r) — (p(r)Q(r). This is close to the expectations fol- 
lowing from the effective-interaction approach Q. Due 
to the boundary condition U(k) — U'(k) cx k 4 at small 
k the mass of the quasiparticles coincides with that of 
the primordial bosons. However, as it has been men- 
tioned in Ref. || , Eq. (|5l]) is the simplest of the possible 
approximations which are fixed by the necessary condi- 
tion that the Bethe-Goldstone equation for ippO 1 ) is re- 
duced to the equation for ip(k) in the limit p — ► 0. These 
approximations lead to the same low-density expansions 
for thermodynamic quantities (the energy, the chemical 
potential, the condensate depletion) in the leading and 
next-to-leading terms. However, difference in the ther- 
modynamics has to appear at small but finite densities 
as well as in the higher-order terms in the density expan- 
sions. The same goes for the microscopical picture given 
by Eqs. (5^) and (|5|). Here the different approximations 
of the relationship between the momentum distribution 
and pair wave functions can lead to some different de- 
tails. In particular, we can not a priori exclude that 
there is a more relevant variant of Eq. ( [5~ij ) which leads 
to U'(k) - U(k) cx k 2 (k 0). In this situation the 
quasiparticle mass would be m* = m/(l + (3n), where 
(3 = \im k ^ Q [U'(k)-U(k)}/T k . 



with T k = T k + n[U'(k) - U(k)}. As to the supracon- 
densate-condensate states, the goal of this paper makes 
it possible not to go into details concerning ip p (r). It 
is sufficient only to use Eq. (55). In the zero-density 
limit n — Eq. (pq) is reduced to the equation ip(k) = 
^°){k) = -U^(k)/(2T k ), which can be rewritten in the 
form of Eq. (|J). So, at small densities all thermodynamic 
quantities can be expressed in terms of the vacuum (out- 
medium) scattering amplitude (k) given by 



U m (k) = / d 3 r</> (0) (r)$(r)exp(-ik-r). 



(59) 



Below it is shown that this feature of Eqs. ( |57j ) and d58| ) 
results directly in Eqs. j65|), j69| ) and (FG). So, the low- 
density energy expansion given by Eq. d7fl) is determined 
in our model beyond any additional assumptions like 
the Landau postulate in the pseudopotential approach. 
Equations ( p7\ ) and (58) yield n k oc l/k and tp(k) cx 1/k 
at small boson momenta. This is totally consistent with 
the well-known "1/& 2 " Bogoliubov theorem [||. It is in- 
teresting that the correct low-momentum behaviour of 
n k and comes from the relation U'(k) — U{k) cx k A 
which follows from Eq. (|54| ) taken in conjunction with 
Eq. (^5|) being a result of the principle of correlation 
weakening. Note that Eqs. (^) and ( p6| ) derived within 
the Bogoliubov model can be obtained from Eqs. ( |57| ) 
and ([38|) by replacing T k and U(k) by T k and <f>(fc), re- 
spectively. So, in what concerns Eqs. (|57j ) and (|58|), the 
situation in our strong-coupling generalization of the Bo- 
goliubov model does look as if we operated with a weakly 



VII. SHORT-RANGE BOSON SPATIAL 
CORRELATIONS 



Now, to elaborate on the picture of the short-range 
boson correlations, let us investigate how the correlation 
hole stipulated by the repulsion between bosons at small 
separations changes due to the influence of the surround- 
ing bosons. At n — > this hole is completely specified 
by the condensate-condensate pair wave function ip(r). 
Exploring how U(k) is expressed in terms of U^°'(k), 
makes it possible to know how tp(r) differs from ip^ (r) at 
small boson separations. Note that the relation connect- 
ing U(k) with [/(°) (k) has been published in our previous 
paper |3| without supporting calculations for reasons of 
space. Let us give these important calculations here. Us- 
ing the definition of U(k) and Eq. (l5S|), for the scattering 
amplitude one can find 



U(k) = - I 



d\ mk-q\)U(q) 
^ ■ If* + 2nf q U{q) 



(60) 



which can be called the in-medium Lippmann-Schwinger 
equation. Let us rewrite Eq. (|60| ) in the form 



U(k) = 



d 3 g $(|k-q|)t/(g) 
(2tt)3 T q 
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where for / we have 



d 3 g / <D(|k-q|)£/(g) _ $(|k-q|)i7( g ) \ 



Performing the "scaling" substitution ( p9| ) in the integral 
and, then, taking the zero-density limit in the integrand, 
for n — > we find 



v? 



■U 3/2 (0). 



(61) 



From Eqs. 



and (61) it now follows that 



U{k)-U {0) {k) =a$(k) 



d 3 q $(|k-q|) 



(2tt) 3 
x[U(q) 



2T„ 



U (0) (?)], (62) 



where U^°\k) obeys Eq. (60) with n = 0, i.e. the 
standard Lippmann-Schwinger equation. Introducing 
the new quantity £(g) = -[[/(g) - f/(°)(g)]/(2Tg), for 
its Fourier transform £(r) we find the equation which 
is nothing else but the Schrodinger equation (^|) with 
ip(°)(r) replaced by a + £(r). As £(r) — > when r — > oo, 
then we can conclude that £(r) = a-!/^ - 1 (r) . Hence, for 
n — * we get 



U(k) ~ C/(°)(fc)(l 



-7(fc,n)^A 

/7T 



(63) 



Here 7(fc,n) — > 1 when n — > 0, and the scattering length 
a is defined by Eq. (|ll]). The derived result for the in- 
medium scattering amplitude U{k) coincides with the 
low-density expansion for the effective potential found 
within the effective-interaction approach at the zero tem- 
perature (see Eq. (4.27) in the review jl8|). This shows 
once more that there are actual parallels between our 
model and approach of Ref. 0. However, these par- 
allels are accompanied by significant differences. First, 
in general the in-medium Lippmann-Schwinger equation 
( |60|) is not a variant of the t-matrix equation being fre- 
quency dependent contrary to Eq. (|60|). Second, Eq. ( |60| ) 
has been found beyond any diagram technique by means 
of the variational procedure whose consequence is that 
the pair wave functions "generating" in-medium scatter- 
ing amplitudes, coincide with the pair wave functions in- 
volved in g{r). On the contrary, it is not true for the 
effective-interaction scheme. It implies the plane waves 
for tfip (r) in the pair distribution function (see the discus- 
sion in Sec. [V|) and certainly goes beyond the plane- wave 
approximation when calculating t-matrix corresponding 
to a pair of particles with a nonzero total momentum. 
Third, in Eqs. (f57|)-(|60"|) we deal with Tk rather than with 
Tfc appearing in the effective-interaction scheme. Fo r 



more information, see also discussion in the paper |10|. 
Now, returning to Eq. (|63|), we can conclude that for 
r ;$ R (R is the radius of the interaction potential $(r), 
for strongly singular potentials R is of the order of the 



scattering length a) and n - 
in-medium renormalization: 



we obtain the following 



<p(r) ~¥> (0) (r)(l + 



(64) 



Thus, the correlation hole coming from the repulsion of 
bosons at small particle separations gets less marked with 
an increase of the density of the surrounding bosons, 
which is mainly the result of the Bose-Einstein statistics. 
For the pair distribution function at small boson densities 
wehaveg(r) cx [(p(°\r)} 2 [see the expression (|S3"1)1 . So, for 
the strongly singular potentials, when ip(°'(r = 0) = 0, 
the correct strong-coupling result g(r = 0) = takes 
place for a dilute Bose gas if Eqs. (g|), m) and ( |5l| ) are 
taken as the basic relations. 

It is interesting to note that for the effective-interaction 
approach 0] the upper cutoff in the momentum space at 
p c ~ H/a is usually made. Thus, one could expect that 
owing to the uncertainty relation the in-medium renor- 
malization is essential when a < r in the real space. 
Equation (|^) shows that this is not the case. As it is 
seen, ip (r) at small separations is really a solution of the 
"bare" Schrodinger equation (see Sec. |) but differs from 

Va/tt)- We also re- 



ip("'(r) by the multiplier (1 
mark that the peculiar overscreening takes place for tp(r) 
when r — > oo. Indeed, Eq. (58) yields 



lim kib(k) 

fc->0 



mU(0) 



The last relation implies that for the condensate- 
condensate pair wave function we have ip( r ) = fi 7 ") — 1 oc 
1/r 2 for r — > oo, in contrast to the "bare" wave function 



^°)(r) 



(0) 



(r) — 1 oc 1/r for r — > oo [see Eq. (M)]- 



VIII. LOW-DENSITY EXPANSIONS 



Now, to verify that a subtle balance of the terms com- 
ing from the short-range particle correlations plays a sig- 
nificant role in the problem of the strong-coupling Bose 
gas, let us calculate low-density expansions of the basic 
thermodynamic quantities. The relation for the conden- 
sate depletion 



n 



d q n q _ 
(2?r) 3 ~n ~~ 30F 



(65) 



can be obtained from Eq. ( J57| ) with the "scaling" substi- 
tution given by Eq. (|39|). 

The low-density expansion for the energy can be de- 
rived in four different ways. 
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A. The chemical potential 

The first way of obtaining the energy expansion deals 
with the chemical potential u and starts from the follow- 
ing relation for a valid in presence of the Bose conden- 
sate §1: 



1 



M = 



dV $(|r - r'|)(V' t (r')V'(r')V'(r))- (66) 



Here ip ( r ) an d V'(r) stand for the Bose field operators. 
This relation follows from the well-known expression for 
an infinitesimal change of the grand canonical potential 
5£l = (S(H — fJ,N)) and the necessary condition of the 
minimum of Q with respect to the order parameter Nq: 
dil(No, fi, T)/dNo = 0, the Hamiltonian depending on 
the number of the condensed particles owing to the sub- 
stitution = ao = \/Nq. Equations (||) and (66) lead 
to M 



n U(0) + V2 



9 £/ q/2 (q/2), 



(2tt) 3 

where U (0) is defined by Eq. ([53]) , and 

U p (k) = J d 3 rv?p(r)$(r)exp(-ik' 



(67) 



(68) 



Using the substitution (|39|) in the integral and taking into 
consideration Eqs. (|5|), ( |63|) and (p5|), we can rewrite 
Eq. @ for n->0as 



, . , . AirTi an , 
(j, = nU(0)(l + C + -") = (1 



(69) 



This, together with the thermodynamic relation u = 
d(ne(n)) /dn, yields the following result: 



150F 



(70) 



known since the familiar paper by Lee and Yang |25| ] and 
found with the binary collision expansion method. 



B. Direct calculation of the energy 

The way of this subsection is direct and starts from the 
expression (||). Inserting Eq. (|^) into Eq. (^) and using 
the substitution (|39| ) in the integral, we can rewrite the 
pair distribution function for n — > in the form 



g(r)~ (l + 2CV 2 (r), 



(71) 



where the relation ( pq) is implied. Note that this expres- 
sion is not valid at sufficiently large r as the boundary 
condition g(r) — > 1 for r — > oo is not satisfied. However, 



here we are not interested in the long-range behaviour of 
g(r) because we use Eq. ( |7l| ) when integrating g(r) mul- 
tiplied by the short-range potential $(r). Equation ( |7l| ) 
makes it possible to represent Eq. dq) for n — ► as 



£ ~-(l + 2C)Z7(0) 



(2tt) 3 



'-(l + 20U(q)^(q) 



(72) 



Taking the term proportional to £ in the integral in 
Eq. (f72l) , we can rewrite it for n — > in the form 



(27T) 



(2tt) s 



U^(q)^°\q) 



-71^- 



6, 



(73) 



where the vacuum scattering amplitude U^~Hq) and the 
characteristic length b are given by Eqs. ([39|) and (|2^), 
respectively. Using the substitution ( |39| ) in the residual 
part of the integral in Eq. (f72|), 



d 3 q 



and taking into account Eqs. (p7|), (p8[), (|63| ) and (|6J 
we arrive at Eq. J70| ) but with the second term multiplied 
by the factor A = 1 — 56/ (8a). Now the question arises 
which variant we should prefer, Eq. @ or Eq. @ with 
the factor A = 1 — 56/(8a), and what is the reason for 
this ambiguous situation? 

To answer this question, let us reconsider the proce- 
dure of calculating e given in this subsection. As it has 
been mentioned earlier, Eqs. (57) and ( p8[ ) used in our 
calculations are consistent with Eq. ( p6| ) rather than with 
Eq. (|5l|). Being characteristic of the Bogoliubov model, 
Eq. (|26[ ) is accurate to the leading order in (n—no)/n and 
differs from Eq. ( pl| ) by the supracondensate-condensate 
term neglected in the Bogoliubov relation. The prob- 
lem of the A-factor turned out to be directly related to 
this term. It can be taken into account by representing 
Eq. (pl|) for n — > in the form 



n k (n k + l) = (l + 2C)V> 2 (fc). 



(74) 



Solving this equation with respect to n k and noticing 
Eq. (|58|), one can obtain 




T k + nU(k)) +2(n 2 U 2 (k) 



T| + 2nT k U{k) 



(75) 



Now, restarting from E q. ( 1 72 ) and making use of the sys- 
tem of Eqs. (J7g) and (|58|) instead of that of Eqs. (57) 
and (|58|), we arrive at Eq. (f70|). The term given by (73) 
is now cancelled due to the correcting term 2(n 2 U 2 (k) 
involved in Eq. (|75|). So, we face a rather c omp licated 
situation: namely, to get the correct result (70) in the 
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direct calculations starting from Eq. (||) , we have to 
abandon Eq. ( p7| ) in favour of Eq. ( J75| ) while for if>(k) we 
can exploit Eq. (|5q). The most important point here is 
the uniform convergence of the integral in Eq. ( f72| ) pro- 
vided Eq. ( f75| ) is used. This allows for employing Eq. ( |75| ) 
together with Eq. ( p8[ ) in spite of the fact that the lat- 
ter has been found in the leading order in (n — no)/n. 
The higher-order corrections to Eq. ([58]) does not influ- 
ence the result of calculating the integral in Eq. (|7|) if 
we limit ourselves to the leading and next-to-leading or- 
ders in na 3 . It is worth noting that replacing Eq. ( |75| ) 
by Eq. ( |57| ) does not influence on Eqs. (|65|) and (^M). So, 
the preliminary result for e found in the paper [M] and 
corresponding to Eq. ( |70| ) with the second term multi- 
plied by the factor A = 1 — 56/ (8a), has to be abandoned 
m favour of Eq. @. 

The analysis fulfilled in this section demonstrates the 
crucial role of the subtle balance of the terms coming 
from the boson scattering (or, in other words, from the 
short-range boson correlations). Disturbance of this fine 
interplay which seems to be insignificant, can never- 
theless lead to wrong conclusions. We stress that the 
strong-coupling model of Ref. || is balanced because it 
takes into consideration the supracondensate-condensate 
scattering waves in both the pair distribution function 
and relation connecting the momentum distribution with 
the pair wave functions. On the contrary, the effective- 
interaction approach is not balanced with respect to the 
supracondensate-condensate scattering waves which are 
missed in the pair distribution function but make a con- 
tribution to the "dressed" potential (see Sec. |y|). Exactly 
this problem is the basis for the so-called ultraviolet di- 
vergence occurring in the effective-interaction approach. 



C. The energy expansion through the 
Hellmann-Feynman theorem 

As it is shown in Sec. [vj the effective-interaction ap- 
proach results in the irrelevant picture of the short-range 
boson correlations. This is why it can not yield the cor- 
rect individual values of the interaction £i„ t and kinetic 
Skin energies by the direct calculations based on Eq. (||). 
Recall that we have the following definitions: these ener- 
gies are defined by 

E ™* = ^(E|^- r il)) = i /d 3 r 7 $(r) 5 (r), (76) 
£km = lv(~^2^ J ~} = I 7Kzk T k-T' ( 77 ) 



2m 



Our model provides the correct short-range behaviour 
of the pair distribution function g{r). So, we can first 
evaluate £mt(7), and, then, obtain the total energy 
(||) by means of the well-known expression often called 
the Hellmann-Feynman theorem which is just the varia- 
tional theorem for the ground state obeying the ./V-body 
Schrodinger equation: 



SE = (6H). 



(79) 



In Eq. (79) 5E and SH are infinitesimal changes of the 
average energy (E = (H)) and the Hamiltonian 



H 



x h 2 v1 1 x - ... 



(80) 



respectively. The relations @, (0), © and © leads 
to the important equations: 



de_ 
dm 



&int — T'o 5 &kin 



From the first expression in Eq. (JSlJ) it follows that 
Jo 7 



(81) 



(82) 



To evaluate £i n t{l) in the leading and next-to-leading 
orders in na 3 , it is convenient to rewrite Eq. (fry) as 



g(r) ~w(na 3 )[^ (0) (r)] 



(83) 



where Eqs. (|64|) and (^35|) are taken into account and 
w(na 3 ) is given by 



w{na 3 ) = 1 + 



(84) 



The range of particle separations for which Eq. (|83j) is 
correct, coincides with that of Eq. (|7l| ) [see the discus- 
sion there]. Keeping in mind Eqs. @, (§|), ©, @ 
and (H^), we obtain 



2Tth 2 n f 1 da , 

/ d 7 — w(na (7)) 

m Jo &y 



2irh 2 n 



da' w(na' ). 



■m jo 

Thus, using Eq. (|84|), we arrive at Eq. (fro| ) again. 

D. The energy expansion through the virial theorem 



where (• • •) stands for the statistical average with respect 
to the ground state, and the auxiliary parameter 7 is the 
coupling constant. The total energy per particle (||) is 
given by the sum of £i nt and e^n at 7 = 1: 



E/N 



(78) 



This method was proposed by Bogoliubov in his orig- 
inal paper || in order to obtain the leading-order term 
in the energy expansion. Here we conside r this method 
in a more general form. As in Sec. VIII C , we start from 
the expression for the pair distribution function g(r, n) 
which is assumed to be a known function of the density 
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n. The basic idea is to derive the differential equation 
for e(n). 

On the one hand, from the virial theorem we get the 
following expression for the pressure 



2 , N n 2 f ■ 
P = -^ekin(n)n — / d" 



d$(r) 



/•— rg(r, n), 



(85) 



where £fcm(n) is given by Eq. (|77|). On the other hand, 
we have the thermodynamic relation 



P = n 



,de(n) 
dn 



(86) 



valid at the zero temperature. Here e(n) is the energy 
per particle (H) which can be written as 



e(n) = e kin {n) 



d 3 r g{r, n)$(r). 



(87) 



The system of Eqs. (85 )-(p7|) yields the differential equa- 
tion for s(n) whose general solution is of the form: 



6 



d 3 r 



d$(r) 
dr 



r + 2$(r) x(r,ra), 



where the function x(r, n) stands for 



X(r, n) = n 



= «a/3 



dn 



,/2/3 



(89) 



and Co is the integration constant. Note that Eqs. ( |88| ) 
and (89) are valid for both the Bose and Fermi systems 
because we have not used the type of the statistics when 
deriving these equations. For a Fermi system the con- 
stant Co is not equal to zero while for a Bose one we 
should put Co = provided the Bosc-Einstein condensa- 
tion takes place. Substituting Eq. fl83| ) in Eq. (|S9|), from 
Eq. (||) we get 



J- 



,2/3 



,w(n'a ) 

1 „>V3 



(90) 



where 
J = 



d 3 r 



d$(r) 
dr 



r + 2$(r) [^ (0) W] 2 = ^ (0) (0)- 



The last equation can be derived after small algebra with 
the help of Eqs. (§)-©. As it is seen, Eq. @, together 
with Eq. leads to the low-density expansion (|70|). 

Thus, all the four ways of calculating e within the 
strong-coupling model developed in Ref. || leads to 
Eq. A 



E. The interaction and kinetic energies 

For any physical quantity there usually exist various 
calculating procedures leading to the same result pro- 
vided the model considered is consistent. By contrast, 



in presence of a thermodynamic inconsistency different 
ways of calculating any thermodynamic quantity are able 
to produce different results ]26|;^] , during which it is of- 
ten happened that one of them is reasonable but others 
are completely inadequate. Such a situation is realized 
when evaluating the interaction and kinetic energies via 
the effective-interaction method. This is demonstrated 
below. 

The interaction 



[6|) and kinetic (77) energies of a 
dilute Bose gas can be evaluated on the basis of the 
Hellmann-Feynman theorem with the help of Eq. (pZC|) . 
Representing this expansion for e in the form 



2nh an 



f(na 3 ) , 



(91) 



and keeping in mind Eqs. wSW) and (p3[), one can derive 



2irTi 2 (a — b)n ( 3 3 d/(na 3 ) 



f{na A ) + ina 6 



m \ d(na 3 ) 

2irh 2 bn ' ' Af(^^ 



(92) 



m 



together with Eqs. (B2J) and <M) yields 



According to Eq. Q f(x) + 128^/(15^), which 

(94) 
(95) 



2irh 2 (a-b)n , 
e int = — I 1 



3^ 



£ hi 



2irh 2 bn 



64 / 3a 

— ^VJia- 3 1 - 

30r V 5& 



As it is seen from Eqs. (|9l|)-(|93]), the terms involving b are 
present in the expressions for the kinetic and interaction 
energy and mutually cancelled in the total energy e. We 
emphasize that the reasoning of this paragraph can be 
fulfilled for both the effective-interaction approach and 
model developed by the present authors. 

Our approach is fully consistent, which makes it possi- 
ble to derive Eqs. (Q) and (^5|) in another way using the 
direct calculations. Indeed, Eq. (|76|) taken at 7 = 1 in 
conjunction with Eqs. ( f2l| ) and (^3[) results in Eq. (|9^). 
Notice that the supracondensate-condensate scattering 
waves make a significant contribution to the next-to- 
leading term of the low-density expansion for £, n t. It 
is not also difficult to find the low-density expansion for 
the kinetic energy ( |95| ) with the help of Eqs. (|||), ( |65| ) 
and @. 

On the contrary, due to the thermodynamic incon- 
sistency the effective-interaction scheme does not allow 
for obtaining Eqs. (Q) and ( |95| ) directly, beyond the 
Hellmann-Feynman theorem taken with Eq. (]7C|). Let 
us show this for the pseudopotential approach discussed 
in Sec. 0. To evaluate the kinetic energy in this case, one 
can start from Eq. ( |77j ) and use the Bogoliubov formula 
( |35|) with the pseudopotential substitution (f47|). In so 
doing the divergent integral J d 3 k/k 2 should be ignored 
[see the discussion in Sec. Similarly, the interaction 
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energy at 7 = 1 can be derived from Eqs. (|25|), (|3£ 
( p6| ) and (|47]) using the same regularization. However, to 
simplify the calculations, we adopt the other way leading 
to the same results and based on the low-density expan- 
sion ( pL^ ) found within the Bogoliubov model. The origi- 
nal Bogoliubov scheme is fully consistent, which implies 
equivalence of different ways of calculating any thermo- 
dynamic quantity. Therefore, we can first find the kinetic 
and interaction energy by using the Hcllmann-Feynman 
theorem together with Eq. ( fi"2|) and, then, replace clq by 
a and a\ by in the derived expressions. From the defi- 
nition ([l6]) it follows that 

da da dai dai 

-f——=m——=a , 7— — = m-— - = 2a x . 

oj am oj dm 

Hence, within the Bogoliubov model we can arrive at 



~ hi- 




(96) 
(97) 



provided Eqs. © and Eq. @ are taken into consid- 
eration. Now, replacing ao by a and substituting for 
di (the latter allows for escaping the ultraviolet diver- 
gence, see Sec. 0), we obtain the following expressions: 



2nh an 



Ski 



1 



m \ 3-s/tt 
2-Kh 2 an 64 , — 



m 



50? 



(98) 
(99) 



which should be compared with the correct results given 
by Eqs. (j94| ) and (P5| ). As it is seen, the sum of the r.h.s. 
of Eqs. (|98j) and (|99D gives the r.h.s. of Eq. ( [70| ) but at the 
expense of the negative value of the kinetic energy (pg|). 
Notice that the Bogoliubov model is free from this non- 
physical feature because ai < 0, which leads to £hm > 
at sufficiently small densities. Thus, Eqs. ( p8| ) and ( |99| ) 
found within the pseudopotential model are inadequate. 
The reason is obvious: the pseudopotential scheme allows 
for restoring the functional dependence on the scattering 
length a in Eqs. (|9^) and (Q) while completely ignores 
the additional length b that cannot be involved in the 
pseudopotential model due to the ultraviolet divergence. 
_Note that in the case of the hard-sphere interaction 
we get a = b from the solution of the Schrodinger 
equation @. Then, in the general case of Eq. (|9|) the 
relations e int = 0, e^ in = e come from Eqs. ( p2[ ) and (p^). 
Thus, we arrive at the interesting property of the Bosc 
gas with the pairwise potential (|l^): namely, although 
the Bose gas is strongly interacting, the interaction en- 
ergy is equal to zero. Hence, the total energy of a dilute 
Bose gas made of the hard spheres is exactly equal to 
the kinetic one. One can see that this result is rather 
general: for the hard-sphere potential ( fil| j the interac- 
tion energy is equal to zero for any density. Indeed, $(r) 



given by ([18] ) can be thought of as the limiting case of 
the repulsive potential 



$(r) 



Vq, r < a, 
0, r > a. 



It is clear that the saturation takes place when Vb 3> £■ 
further increase of the parameter Vb does not change the 
energy per particle e. Hence, according to Eqs. (76|) and 
(|si|), ei n t = because de/dj = at 7 = 1 in the limit 
Vb — > +00. Notice that even taken in the order linear 
in the density n, Eqs. ( p8| ) and ( |99| ) lead to the opposite 
case £mt — s, £kin — 0. This incorrect redistribution 
of the energy of a dilute Bose gas in the pseudopoten- 
tial approach is also remarked in the paper |2g| l where 
the leading order of the low-density energy expansion is 
considered. 

Note that the relation Eq. (|2^) enables us to obtain 
the length b in an experimental way, by the isotopic shift 
of the scattering length a: 



a(l 



d In a 
<91nm 



a 1 



Aa m 



Hence, we are able to evaluate the interaction (J9J) and 
kinetic (^5|) energies per particle via the quantities which 
can be found experimentally. 



IX. CONCLUSION 

In conclusion, we remark that this paper concerns the 
thermodynamics of a dilute Bose gas with a strongly re- 
pulsive interaction in the leading and next-to-leading or- 
ders of the low-density expansion. The strong-coupling 
generalization of the Bogoliubov model proposed by the 
present authors is shown to reproduce the result ( |70| ) 
of Lee and Yang |2^] found via the binary collision ex- 
pansion method. Contrary to the effective-interaction 
approach of Ref. Q, the model considered in this pa- 
per is thermodynamically consistent and free of the ul- 
traviolet divergences. These advantages are due to ac- 
curate treatment of the short-range spatial boson corre- 
lations whose picture is inadequate within the effective- 
interaction scheme. The present paper thus demonstrates 
that the effective-interaction scheme which is reduced to 
the Bogoliubov model with the effective pairwise poten- 
tial, is not acceptable for investigating a dilute strongly 
interacting Bose gas. In addition to the arguments men- 
tioned above, this also follows from the results for the 
kinetic and interaction energies first found in this paper. 

In some sense the strong-coupling model discussed can 
be considered as a generalization of the Brueckner ap- 
proach taken in its representation given by Bethe and 
Goldstone [^9). The new essential point is that the in- 
medium pair wave functions are calculated in conjunction 
with the particle momentum distribution on the basis of 
the variational procedure. So, to go further, additional 
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investigations should be fulfilled to establish more accu- 
rate approximations of the relation connecting the boson 
momentum distribution with the scattering parts of the 
in-medium pair wave functions. In particular, this im- 
provement is needed to clarify to what extent the correct 
spectrum of elementary excitations in a dilute Bose gas 
differs from the well-known prediction of the effective- 
interaction approach. Of course, the region of interme- 
diate momenta is implied rather than the linear phonon 
sector which should be the same according to the thermo- 
dynamic prescription. This problem is tightly related to 
investigation of the long-range spatial boson correlations 
beyond the effective-interaction approach. 
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